The quasielastic ^H(e, e'p)n reaction at high recoil momenta 
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The ^H(e, e'p)n cross section was measured in Hall A of the Thomas Jefferson National Accelerator 
Facility (JLab) in quasielastic kinematics {x = 0.96) at a four-momentum transfer squared, Q^=0.67 
(GeV/c)'^. The experiment was performed in fixed electron kinematics for recoil momenta from zero 
to 550 MeV/c. Though the measured cross section deviates by 1 — 2(t from a state-of-the-art 
calculation at low recoil momenta, it agrees at high recoil momenta where final state interactions 
(FSI) are predicted to be large. 
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The deuteron, as the only bound two-nucleon system, 
represents the simplest manifestation of the nuclear force 
and is therefore the natural starting point for investiga- 
tion of the nature of the nuclear electromagnetic current. 
The applicability of reaction models for complex nuclei 
can be gauged by the success of these models in repro- 
ducing scattering observables on the deuteron. Further, 
by studying the deuteron in extreme kinematics where its 
short range structure is emphasized, one may determine 
to what extent its description in terms of nucleon/meson 
degrees of freedom must be supplemented by inclusion of 
explicit quark/gluon effects. This issue is of fundamental 
importance to nuclear physics. Finally, understanding 
the deuteron is required for the interpretation of a va- 
riety of experiments using the deuteron as an effective 
neutron target. 

Coincidence ^H(e, e'p)n reactions are particularly well 
suited to nucleon-nucleon (JVJV) interaction studies. Be- 
low pion threshold, the final state is completely speci- 
fied. Further, by judicious choice of kinematics one can 
emphasize various aspects of the reaction dynamics [|]. 
Accessing the deuteron's short range structure implies 
measurements at high recoil momentum, pr (i.e. the mo- 
mentum of the undetected recoiling neutron). However, 
for certain kinematics, final state interactions (FSI) can 
change the cross section by an order of magnitude or 
more at high pr (see, e.g.. Fig. |l| below). These large 



effects result from strength at low initial proton momen- 
tum feeding high pr due to np rescattering in the final 
state. Models which succeed under such stringent tests 
may then be applied with some confidence to infer as- 
pects of the deuteron short distance structure in kine- 
matics where Pr is large but FSI are minimized. Such a 
situation is expected at large x {x ~ Q'^/2mLu) in "par- 
allel kinematics" , where the proton is detected along the 
direction of ||] . 

Although there is a substantial body of data on this 
reaction from facilities other than JLab, including unsep- 
arated cross sections [3, I , p|, pi a s well as separations of 
response functions 0, 3, 3, 0,^, [l3[ |l|| , various 
limitations made it difficult to disentangle aspects of the 
reaction mechanism from the deuteron short range struc- 
ture. For the cross section measurements, limitations 
in the accelerator energies available at Bates, Saclay, 
NIKHEF and Mainz frustrated such attempts by forc- 
ing measurement of very high recoil momenta to energy 
transfers far above the quasielastic peak. Thus, for the 
Turck-Chieze and highp^ Mainz data, the kinemat- 
ics were in the A-region where lack of knowledge of the 
reaction mechanism made it difficult to deduce aspects 
of the deuteron structure. Although the Mainz measure- 
ment sampled Pr up to 928 MeV/c, the kinematics actu- 
ally imply that the bulk of the cross section arises from 
interaction with the neutron, leaving the detected pro- 
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ton as a spectator. Further, since the kinematics were in 
the A- region of the inclusive (e, e') spectrum, the inclu- 
sion of virtual nucleon excitations was required to obtain 
agreement with the data. Although the energy limita- 
tion is not shared by SLAC, the maximum current and 
duty factor restricted the range of recoil momenta there 
as well. In contrast, JLab's combination of high beam 
energy, current and duty factor allows examining large 
recoil momenta at or even below quasielastic kinemat- 
ics (i.e. X > 1), making the extraction of the deuteron 
structure less model-dependent. 

The ^H(e, e'p)n separation experiments have also been 
restricted in kinematical coverage for the reasons stated 
above. Nonetheless, they have revealed gaps in our un- 
derstanding. Various calculations have difhculty repro- 
ducing both Rl and Rt H, |11|, E3|. The Rlt response 
and related asymmetry ^ |l^ |ll[ |l^ indicate the 
need for relativistic treatments but problems still exist 
in reproducing the data. For momentum transfers up 
to 1 (GeV/c)^ a nonrelativistic calculation supple- 

mented with the most important leading order relativis- 
tic contributions has been shown to agree quite well with 
a covariant approach |l6j. More recent calculations of 
Jeschonnek, Donnelly and Van Orden Q suggest that 
the bulk of the relativistic effects may be in the nucleon 
current operator, as opposed to the nuclear dynamics. 
They obtained good agreement between a manifestly co- 
variant calculation involving the Gross equation and 
a calculation using a nonrelativistic wave function and a 
relativistic nucleon current operator. In contrast, a cal- 
culation with a nonrelativistic wave function and current 
operator drastically failed to reproduce the results of the 
covariant calculation. These calculations highlight the 
importance of incorporating relativity properly. Also, if 
the relativistic effects are mostly in the current opera- 
tor, they can be incorporated for heavier nuclei using the 
same approach. Due to its more tractable nature, the 
deuteron provides a testing ground for such an approach. 

Our experiment (Experiment E94-004) consisted of 
measuring the unseparated cross section in quasielastic 
kinematics out to high pr- It was performed in Hall 
A of JLab using the high resolution spectrometer pair. 
A detailed description of the Hall A instrumentation is 
currently being prepared [p^ . Electrons from the accel- 
erator were incident on a 15 cm long liquid deuterium 
target. Scattered electrons and knockout protons were 
detected in the "left" (formerly "electron" ) spectrometer 
and "right" (formerly "hadron") spectrometer, respec- 
tively, each with nominally 6.0 msr collimators at their 
front ends. Each spectrometer was equipped with its 
standard detector package consisting of a pair of VDC's 
for track reconstruction and a scintillator array for trig- 
ger definition. In addition, the electron spectrometer in- 
cluded an atmospheric pressure CO2 threshold Cerenkov 
detector for 7r~/e discrimination. 

Kinematics are given in Table | The beam energy 



TABLE I: Kinematics (central values). The electron kine- 
matics were fixed throughout the experiment with incident 
and scattered energies of £ = 3.110 GeV and E' = 2.742 
GeV, respectively, and scattering angle of 6^ = 16.06°, lead- 
ing to central values of the four-momentum transfer squared 
of = 0.665 (GeV/c)^ and x = j^muj = 0.964 (where 
is the electron energy transfer and m is the proton mass). 
Listed below are the recoil momentum (pr), the proton mo- 
mentum (p), and the proton angle (Qy) relative to the beam 
direction. All quantities are expressed in the laboratory sys- 
tem. Also shown are the average luminosities corrected for 
estimated beam related target density changes. 
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300 
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500 
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-92.78 


3.67E-f38 



and electron spectrometer angle and central momentum 
were kept fixed throughout the experiment at values cor- 
responding to the top of the quasielastic peak. The recoil 
momentum was varied by changing the proton spectrom- 
eter angle starting from parallel kinematics (i.e. protons 
detected along the three-momentum transfer, q) where 

= to more backward angles resulting in a maximum 
central value of = 500 MeV/c. The proton spectrome- 
ter momentum setting was varied in conjunction with the 
angle so that the central kinematics were fixed near the 
deuteron breakup energy of 2.2 MeV. Obstruction from 
one of the scattering chamber support posts resulted in 
a fairly wide spacing between the last two settings and a 
gap in the measured spectrum vs. Pr- 

Measurement of the elastic ^H(e,e'p) reaction, taken 
with the spectrometers at the ^H(e, e'p)n pr = Q setting, 
served as a normalization check. The measured yield was 
compared to a simulation |po| using the Simon et al. ]2l|] 
parameterization for Gm^ and the Gep/Gm^ ratio mea- 
sured by Jones et al. H^. The simulation included accep- 
tance averaging as well as radiative folding |23| . In order 
to remove the ill-defined acceptance edges a technique 
involving "i?- functions" was used |2j]. Through the R- 
functions, a multi-dimensional contour in the space of 
the target variables and equidistant from a pre-defined 
boundary was defined and all events outside the contour 
were rejected. The same contour was used in the simu- 
lations and for the "^H(e, e'p) and ^H(e, e'p)n data. The 
ratio of integrated yield for the simulation to that for the 
^H(e, e'p) data was 1.054, amounting to a 5.4% correction 
for the ^H(e, e'p)n data. 

For the ^H(e, e'p)n data, beam currents ranged from 
10 /iA to 100 /iA. Since the electron kinematics were 
fixed, the electron arm served as a measure of the prod- 



10" 



uct of electronic deadtime and target density; the elec- 
tronic deadtime for the hadron arm was assumed to be 
negligible, since its trigger rate never exceeded 10 kHz. 
Computer deadtime was determined from the ratio of 
coincidence raw triggers to recorded events. In order to 
improve the reals-to-accidentals ratio, especially impor- 
tant at high pr, a consistent vertex from both spectrom- 
eters was required and a cut on the missing mass was im- 
posed. Finally, a cut on the summed analog signal from 
the Cerenkov detector was used to reduce the already 
small contamination from 7r~ events in the electron arm. 

The aperture/magnetic model of each spectrometer 
was tested by measuring a series of "white" spectra, 
scanned in overlapping momentum steps. The relative 
spectrometer acceptance was then extracted by an itera- 
tive procedure. For the i?-function cuts described above, 
the results were in excellent agreement with the simu- 
lated phase space. 

The total systematic uncertainty in the cross sections 
was estimated to be roughly 8%, nearly independent of 
Pr, except near zero where it grows sharply due to the 
shrinking of the phase space volume. The kinematic 
related uncertainties were estimated by computing the 
cross section at the center of a given bin in p^, further 
weighted by p^ to roughly account for the phase space 
volume, and then making variations of each kinematic 
quantity in turn. The quantities were varied in a corre- 
lated fashion, constrained by the elastic ^H(e, e'p) kine- 
matics and by independent measurement of the beam en- 
ergy, obtained by measuring the position of the beam at 
a point of high dispersion in an eight dipole arc section. 
The quadrature sum of all kinematics related uncertain- 
ties was less than 2.5% for all recoil momenta not too 
close to zero. Other uncertainties were estimated to be: 
2.0% (electronic deadtime x target density), 2.0% (elec- 
tron arm solid angle; not accounted for in the "'^H(e, e'p) 
normalization where the acceptance was limited by the 
proton arm), 1.0% (beam charge relative to normaliza- 
tion), 3.0% (radiative correction) and 5.8% (normaliza- 
tion, consisting of 3.6% from uncertainty in the elas- 
tic form factors, 2.0% uncertainty from kinematics and 
4.1% from other uncertainties specific to the normaliza- 
tion measurement). Adding these uncertainties quadrat- 
ically yields a value of 7.2% taken to be independent of 

Pr- 

The results for data along with various calculations 
are shown in Fig. ^ The top panel shows the radiatively 
corrected "reduced" cross section: 
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where if is a kinematic factor, crcci is the half-off-shell 
electron proton cross section of de Forest Q and /rec is a 
recoil factor which arises from the integration over miss- 
ing mass. This division removes most of the kinematic 
dependence, except for the Pr dependence, and results 
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FIG. 1: Top panel: The reduced ^H(e,e'p) n cross section for 
this experiment along with various model calculations (see the 
text for details). Bottom panel: Cross sections for data and 
calculations shown as percentage deviations from the "full" 
calculation of Arenhovel. Also shown is the systematic er- 
ror band (ilcr), arbitrarily placed on the vertical axis. This 
error band contains an overall 7.2% contribution added in 
quadrature with the kinematic contribution, the latter vary- 
ing slightly with p,-. 



in a smooth spectrum vs. Pr- The bottom panel shows 
the relative deviation of data and theoretical predictions 
from the "full" calculation of Arenhovel |^ (described 
below). Also shown in the bottom panel is the system- 
atic error band, arbitrarily placed vertically. The data 
were radiatively corrected by multiplying the measured 
cross section in each p,. bin by the ratio of yield without 
and with radiative effects (internal and external), both 
estimated using Arenhovel's full calculation. The cal- 
culation was radiatively folded using the model of Borie 
and Drechsel . Both data and simulations were cut on 
the missing mass from —3.5 MeV to 10.5 MeV, resulting 
in a weak dependence on pr for the radiative correction 
factor. A straight-line fit to this dependence was used. 

All of the models were acceptance averaged using 
MCEEP pOl 



the calculations of Arenhovel were per- 
formed on a grid over the experimental acceptance and 
interpolated for each event, whereas the calculations 
of Jeschonnek were incorporated directly into the 
Monte Carlo simulation program as a subroutine pack- 
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age. The Jeschonnek calculation was in the plane wave 
Born approximation (PWBA) and used a fully relativistic 
single-nucleon current operator with an alternate three- 
pole parameterization of the MMD nucleon form factors 
p8[ and the Argonne V18 two-body interaction ||2^. The 
calculations of Arenhovel included relativistic contribu- 
tions of leading order in p/m to the kinematic wave func- 
tion boost and to the nucleon current. The Bonn r-space 
JVJV potential Q and dipole nucleon form factors were 
used. The various curves are for PWBA, distorted wave 
Born approximation (DWBA, which includes FSI), and 
the "full" calculation which also includes non-nucleonic 
currents: meson exchange currents and virtual nucleonic 
excitations. 

The two PWBA calculations are reasonably close to 
each other, but deviate at high pr, presumably largely 
due to the different JVJV potentials employed. For pr > 
300 MeV/c the PWBA fails completely, whereas the full 
calculation including FSI and significant contributions 
from non-nucleonic degrees of freedom results in satisfac- 
tory agreement with the data. For pr < 100 MeV/c the 
full theory deviates from the data by roughly la growing 
to about 2(7 at Pr ~ 200 — 300 MeV/c, adding statistical 
and systematic errors quadratically. It is important to 
resolve any possible discrepancies in this low pr region, 
especially in light of some of the neutron form factor 
measurements which exploit the deuteron in this kine- 
matic region. Clearly, additional, higher precision mea- 
surements would be helpful in clarifying this situation. 
An experiment to systematically study this reaction over 
a broad kinematical range has already been proposed and 
conditionally approved by the JLab Program Advisory 
Committee This future experiment promises signif- 
icantly smaller systematic errors, based on recently ac- 
quired experience with the experimental apparatus and 
on inclusion of a larger set of experimental cross checks 
and calibrations. 

In summary, we have measured the ^H(e, e'p)n cross 
section at — 0.67 (GeV/c)^ and x — 0.96 for recoil 
momenta up to 550 MeV/c. The data indicate large FSI 
and substantial non-nucleonic effects at high pr . The full 
calculation overestimates the cross section by 1 — 2cr at 
small recoil momenta. 
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